Phase Separation of the Two-Dimensional t — J model 
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The boundary of phase separation of the two-dimensional t — J model is investigated by the power- 
Lanczos method and Maxwell construction. The method is similar to a variational approach and 
it determines the lower bound of the phase separation boundary with J^/t = 0.6 ±0.1 in the limit 
rie ~ 1. In the physical interesting regime of high Tc superconductors where 0.3 < J /t < 0.5 there 
is no phase separation. 



PACS numbers: 74.20.-z, 71.27.-|-a, 74.25.Dw 

It is believed that the main physical properties of the 
high-temperature superconductors can be described by 
the two-dimensional (2D) t— J model on a square lattice. 
The Hamiltonian is: 
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where < i,j > is the nearest-neighbor pairs and ~ 
CiCT(l — 'T^i.-o-)- In this model the two terms compete with 
each other. The kinetic term favors the phase which the 
electrons are homogeneously distributed in the plane to 
minimize the kinetic energy. While the exchange term 
attracts the electrons together to lower the magnetic en- 
ergy. It is easy to see that for very large J jt the system 
will phase separate into a hole-rich region and a region 
without holes to maximize the magnetic energy gain. 

There are experimental evidences as well as theoretical 
studies that indicate phase separation and superconduc- 
tivity are closely related. It is even argued that the driv- 
ing mechanism of superconductivity is the same as that 
of phase separation [l] or superconductivity comes from 
the frustrated phase separation |^ . Hence it is extremely 
important to determine the phase separation boundary 
of the 2D t — J model to resolve these issues. This paper 
reports our findings of the phase separation boundary. 

Experimentally, phase separation of the superconduct- 
ing La2Cu04+5 compound are observed by several mea- 
surements [^-^. The compound phase separates for 
0.01 < (5 < 0.06 below Tp^ w SOOiC into the nearly sto- 
ichiometric antiferromagnetic La2Cu04+5j with 8\ less 
than 0.02 and Neel temperature Tjv ~ 250i4', and a 
metallic superconducting oxygen-rich phase La2Cu04+52 
with (52 ~ 0.06 with ~ MK. The Sr doped compound 
La2-2:Sr2:Cu04+5 also phase separates for x < 0.03 into 
superconducting h8L2-xSrxGu04+s' {S' ~ 0.08) and non- 
superconducting La2-a;Sra;Cu04+5" {S" ~ 0.00) phases 
[0 . Recent muon spin resonance and nuclear quadrupole 
resonance experiments |P-p^ on La2_a;Sra;Cu04 also in- 
dicate that the doped holes were inhomogeneously dis- 



tributed mesoscopically and segregated into walls sepa- 
rating the hole-poor antiferromagnetic domains. 

Theoretically, there are conflicting results. The first 
important paper on this issue is by Emery et al. [0. 
They used the exact diagonalization (ED) to study the 
4x4 cluster. Using Maxwell construction they claimed 
that phase separation occurs for all values of J/t. This 
result is contradictory to the later calculations by using 
quantum Monte Carlo (QMC) ll|] and ED ^ on the 
Hubbard model, which should be consistent with the t — J 
model for small J/t. Putikka et al. studied this problem 
using the high-temperature series expansion and found 
phase separation at T=0 for J/t lying above a line ex- 
tending from J/t = 3.8 at zero filling to J/t = 1.2 at half 
filling |lj] . Prelovsek et al. jl^] calculated the two-point 
and four-point density correlations using ED on clusters 
of size 18 and 20 sites. They found the two-hole bound 
state for J/t > 0.2. For J/t > 1.5 the holes form do- 
main walls along (1,0) or (0,1) direction, and phase sepa- 
rate into a hole-rich and a hole-free phase for even larger 
J/t > 2.5. Hellberg et al. determined very accurately 
that the critical J/t for phase separation at low electron 
density limit is J/t = 3.4367 Poilblanc calculated 

the energy of 2 and 4 holes by ED on several clusters 
up to 26 sites. The phase diagram includes a liquid of 
d-wave hole pairs for J/t > 0.2, a liquid of hole droplets 
(quartets) for larger J/t > 0.5, and at even larger J/t, 
an instability towards phase separation ||lj]. Yokoyama 
et al. investigated the phase diagram by the variational 
Monte Carlo (VMC) method ||l||. The critical J/t for 
phase separation at the high density limit they found is 
1.5, which is consistent with Putikka et at. 

Most recently Hellberg and Manousakis |l^ investi- 
gated this problem by the Green Function Monte Carlo 
(GFMC) method and Maxwell construction for larger 
clusters. Their phase diagram is similar to Emery et al. 
]ll[ . They conclude that the t—J model phase separates 
for all values of J/t in the low doping regime. 

The theoretical results of different groups discussed 
above are consistent at the large J/t and low electron 
density region. But unfortunately, in the interesting 
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physical regime of high Tc superconductors, 0.3 < J/t < 
0.5 and high electron density 0.75 < rig < 0.95, they are 
in disagreement. We have used the power-Lanczos (PL) 
method [^J21| to obtain the best estimate of the ground 
state energy in this physical regime for the largest clus- 
ter (82 sites) that have been studied so far. Based on 
the variational argument we show that there is no phase 
separation in this physical regime |22]. 
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FIG. 1. (a): Typical plots of energy per site vs powers for 
J/t = 0.6, rie = 32/36 (open circles) and = 44/50 (full 
circles) — 74/82 (open triangles); (b) Energy per site as 
a function of electronic density for J/t = 0.6 with different 
cluster sizes. Diamonds are the exact result of 16 sites. Open 
circles are for 36 sites and full circles are for 50 sites, both 
are obtained by Phlpower^e- Triangles are for 82 sites with 

The ground state energy of the Hamiltonian of equa- 
tion (1) is calculated by using the PL method. The PL 
method we used is similar to the GFMC method but 
without using importance sampling and the fixed node 
approximation. The method is essentially a variational 
approach. Applying more powers to a trial wave function 
implies a better approximation of the ground state wave 
function. Details of the method are discussed in Ref. 
pO| . The trial wave functions we used are the optimized 
Gutzwillcr wave functions, resonating valence bond state 
(RVB) [p3| , and RVB with antiferromagnetic long range 
order p4[ . In Fig. 1(a) energy per site is plotted as a 
function of power for J/t = 0.6 and three different den- 
sities: fie = 32/36 (open circles) and Ue = 44/50 (full 
circles) rig = 74/82 (open triangles). Error bar is shown 
only when it is larger than the symbol. We also com- 



pared our energy of J/t = 1 for 50/64 with the result of 
high-temperature series expansion psf . The best energy 
we get is -1.183(2) while the high-temperature expansion 
result is -1.20(2). They are well in agreement. In Fig. 
1(b) we show the best energies we are able to obtain for 
clusters with 36, 50 and 82 sites as a function of elec- 
tronic density. For comparison we also show the exact 
energies of 16 sites ||l^. Energies are little lower for the 
smaller clusters. For 50 and 82 sites, there seems to be 
very little finite size effect. The energy per site is a fairly 
smooth function of density. We do not find large effect 
due to different Fermi surface topology in the physical 
regime. 

To find the phase separation boundary by using 
Maxwell construction we are interested in the variation 
of the slopes in figures like Fig. 1(b). In other words 
we are interested in the second derivative of energy with 
respect to the electronic density, or the inverse compress- 
ibility. It turns out that there is a systematic variation of 
this quantity as the energy approaches the ground state 
or as the power increases in our PL method. Although 
in the physical regime most of our best data have not 
yet converged to the exact ground state, this systematic 
variation is enough for us to determine the lower bound 
of the phase separation boundary. 

It is difficult to read out the slope variation from fig- 
ures like Fig. 1(b), as the curve is almost a straight line for 
He > 0.85. Therefore we shall follow Emery et al. ||ll| by 
examining another quantity. In the one-dimensional t~ J 
model the phase separated state contains an electron-free 
and a electron-rich phases. However, it phase separates 
into a hole- free phase, i.e., the antiferromagnetic Heisen- 
berg island, and a hole-rich phase in the two-dimensional 
t — J model. Thus the energy of the phase separated 
state is in the form: 



E^iNs- N)eH + Nen 



(2) 



where Ng is the total number of sites and N is the number 
of sites in the hole-rich phase, en = 1.169 J denotes the 
Heisenberg energy per site . And eh is energy per site 
in the uniform hole-rich phase, which is a function of the 
hole density in this phase x = N^/N. Nh is the number 
of holes. E can be rearranged into the form: 



where 



E = NsCH + Nhe{x) 



e{x) = [-CH + ehix)]/x 



(3) 



(4) 



If e(x) of a particular J/t has a minimum at a; = Xm and 
the hole density of the total system is smaller than Xm, 
the system will adjust the size of the hole-rich phase N 
such that Xm is equal to N^/N and it minimizes the total 
energy in Eq.(^). Since Ns, en, and Nh are all constants, 
the total energy is minimized as e{x) is minimized. Thus 
Xm is the critical density for phase separation at this J/t. 
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We calculated e{x) from the energy of the uniform 
states eh{x) by the PL method and found the minimum 
of e{x) on 6 X 6, VSO x \/50, and ^/82 x ^82 clusters 
for several densities and J/t. It is very difficult to get 
the converged ground state energy in the physical regime 
due to the sign problem. After we have found the op- 
timized wave function in the VMC calculation we used 
the PL method to project the trial wave function onto 
the ground state systematically. The PL-1 power=4 (for 
82 sites) or PL-1 power=6 (for 50 and 36 sites) energy is 
used here as the eh{x). It is about 2 ~ 4 percent lower 
than the variational energy. We estimate the difference 
between the best PL energy is within one or two percent- 
age of the true ground state energy. 
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FIG. 2. eix) vs hole density x for (a) J/t = 0.4, (b) 
J/t = 0.6 and (c)J/i = 1.5 for several powers: FLO- VMC 
(open circles), PLl-VMC (full circles), PLl-power=2 (open 
triangles), PLl-power=4 (full triangles), and FLl-power=6 
(open squares), (d) J/t = 0.4 for close shells for different size 
of lattices, 74/82, 42/50, and 50/64. 

e{x) vs X ~ 1 — Tie calculated on 50 sites for J/t— OA, 
0.6 and 1.5 is shown in Fig|^(a)-|2|(c), respectively. It is 



interesting to note the trend of the shift of e{x) with 
powers. For J/t = 0.4 (Fig.| (a)), at the VMC level, the 
minimum of e{x) is at x^ — 0.16. It shifts to x = 0.04 
(the minimum hole density we calculated for this cluster) 
immediately after the first order Lanczos improvement 
(PLl-VMC) and stays at the density up to 6 powers. 
For J/t = 0.6 (Fig.|(b)), Xm shifts from x = 0.2 (VMC) 
to X = 0.16 (PLl-VMC) and to a; = 0.08 (PLl-power=6) 
at last. For J/t = 1.5 (Fig.H (c)), Xm shifts from x = 
0.48 (PLl-VMC) to a; = 0.4 (PLl-power=2) and to a; = 
0.36 (PLl-power=6) at last. It is clear that Xm shifts 
monotonically toward a smaller value when the energy 
moves closer to the ground state. 

The results presented in Fig.||(a)-||(c) are calculated 
with a fixed lattice size and different electron numbers. 
Hence Fermi surfaces have different shapes and, in par- 
ticular, there are open and closed shells. It has been 
argued [^9) that comparing energies obtained for these 
different Fermi surfaces might be inaccurate. To exam- 
ine this argument carefully, we have compared systems 
with closed shell Fermi surfaces only. In Fig||(d) e{x) 
calculated from close shells of different size of lattices 
for J/t = 0.4 shows similar behavior as Fig.^(a). The 
minimum of e{x) shifts toward smaller hole density. The 
trend of Xm moving with increasing power is the same 
for both close and open shells. Hence the shell effect is 
not important here. 
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FIG. 3. Phase separation boundary on the phase diagram 
of the two-dimensional t — J model evaluated by ED on the 
4x4 lattice[ll] (open diamonds), by the high-temperature 
series expansion[14] (dashed line) by the GFMC method[19] 
(dotted Une), and the PL method on 36 sites (full triangles), 
50 sites (full circles), and 82 sites(full square). The phase 
boundary determined by the VMC method for 36 sites (open 
triangles) and 50 sites (open circles) are shown in the inset. 

In Fig.^ we show the phase separation boundary deter- 
mined by the best Xm- The PLl-power=6 phase bound- 
aries of 36 sites and 50 sites are shown as full trian- 
gles and full circles, respectively. Also some of the PLl- 
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power=4 data of 82 sites are also shown as full squares. 
For J/t = 0.6 the error bars of the e{x) for Ue = 80/82, 
78/82, and 76/82 are larger than the difference of these 
three e(x), thus error bars of shown in the fieure 

near these electron densities. 

The dashed line in FigJs is the result of high temper- 
ature series expansion p4[]. Similar result is obtained 
by the variational study |[18[. They assumed the sys- 
tem separates into a hole-free Heisenberg antiferromag- 
net and an electron-free vacuum state. This overesti- 
mates the energy required for the phase-separated state, 
since electrons can "evaporate" from the Heisenberg is- 
land to gain energy. Their critical Jc/t « 1.2 is larger 
than our Jc/t « 0.6. Similar argument was also given by 
Hellberg and Manousakis jl^]. 

Our estimate of the Jc/t = 0.6 ±0.1 is actually a lower 
bound. The exact phase separation boundary should be 
to the right of our result in Fig.|^. When we use much 
poorer estimate of the ground state energy as our VMC 
result, the phase boundary is shifted lower. This is shown 
in the inset of FigJ|. The VMC results of 36 sites (open 
triangles) and 50 sites (open circles) show a much smaller 

0.6, 




FIG. 4. e{x) difference between PLl-power=6 and 
PLO-VMC for 50 sites. The values are proportional to the 
area of the circles. 

Another way to understand this argument of lower 
bound is to examine the variation of e{x) with power. In 
Fig.^ we show the change of e{x) between PLl-power=6 
and VMC for 50 sites. The values are proportional to 
the area of the circles. Because of the 1/x factor in Eq. 
(4) , the smaller the hole density the more improvement of 
e(x) will likely occur. Because of the variational nature 
of the PL method, the larger the improvement observed 
between VMC and PLl-power=6 the larger the differ- 
ence between the exact result and PLl-power=6 will be. 
Hence, once the minimum Xm is at the lowest hole density 
such as J/t = 0.4 in FigJ^, better estimate of the ground 
state energy by applying more powers will not change 
the minimum to higher hole density. Based on this argu- 
ment we are confident to conclude that there is no phase 



separation in the physical regime where 0.3 < J/t < 0.5. 

We have found that for J/t < 0.5, the minimum of 
e{x) is always at two holes for clusters of different sizes 
(16, 36, 50 and 82). As argued by Dagotto et al, this 
might indicate a two-hole bound state Q but not phase 
separation. If there were phase separation, the x„i would 
be at the same (or nearby) density rather than the same 
number of holes. 

It is also interesting to note that in Fig.|^(b) , for J/t = 
0.6 the minimum Xm seems to be at 4 holes instead of 2 
holes. This is observed for both 36 and 50 sites. It seems 
to be quite consistent with a recent claim by Poilblanc 
|l7| that there is a phase with quartets for 0.5 < J/t < 
0.8. But our data is not accurate enough for 82 sites to 
make a more definite conclusion. 

Recently Hellberg and Manousakis ||l^ have used 
GFMC to determine the phase separation boundary. The 
phase boundary they reported (dotted line in FigJ|) is 
similar to our variational boundary (see the inset of 
Fig.^). Without knowing details of their calculation we 
cannot completely understand this discrepancy. A pos- 
sible clue is that they might not have obtained lower 
enough energy in the high electron density regime. As 
shown in FigJ|, in particular Fig.^(d), until the energy is 
lower enough to be closer to the ground state, it is very 
easy to to make the conclusion that there is a minimum 
of e{x) at a finite hole density. 

In summary, we determined the phase separation 
boundary by the PL method and Maxwell construction. 
We have studied various size of clusters and densities of 
holes. The largest cluster studied is 2 holes in an 82 site 
lattice. Using the variational nature of the PL method 
and the systematic variation of the energy as a function of 
hole density we conclude that the critical Jc/t for phase 
separation in the low hole density limit is at least « 0.6. 
There is no phase separation in the physical regime. 

It should be pointed out that the result reported above 
are obtained by assuming the hole-rich region in the 
phase separated state has a uniform hole density. We 
have not yet considered more exotic possibilities such as 
the stripe phase ^5 27,2^]. 
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